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The approximate tree decays B~ — > 7r~7r°/p - p may serve as benchmark channels for testing the 
various theoretical descriptions of the strong interaction dynamics in hadronic B meson decays. 
The ratios of hadronic and differential semileptonic B — > nlv/ plv decay rates at maximum recoil 
provide particularly clean probes of the QCD dynamics. We confront the recent NNLO calculation 
in the QCD factorization framework with experimental data and find support for the factorization 
assumption. A detailed analysis of all tree-dominated B —* nn/np/pp decay modes seems to favour 
somewhat enhanced colour-suppressed amplitudes, which may be accommodated in QCD factoriza- 
tion by a small value of the first inverse moment of the B meson light-cone distribution amplitude, 
Xb — 250 MeV. Precise measurements of the B — ► plv spectrum could help to clarify this point. 
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I. INTRODUCTION 

A wealth of observables at current and future B physics 
experiments is related to exclusive hadronic decay modes. 
B decays into a pair of light (charmless) mesons are of 
particular phenomenological interest as they are medi- 
ated by rare flavour-changing b — > q (q = it, d, s) quark 
transitions and the interference of several weak decay am- 
plitudes may induce sizeable CP-violating effects. 

The complicated strong interaction dynamics in ha- 
dronic decays poses a serious challenge for accurate 
theoretical predictions. In recent years systematic meth- 
ods have been developed, which are based on the factor- 
ization of short- and long-distance effects in the heavy 
quark limit m& ^> Aqcd- The theoretical concepts are 
known as QCD factorization (QCDF) soft-collincar 
effective theory (SCET) @ and the pQCD approach Q. 

In this letter we consider the decays B~ — ► ir~ ir° / p~ p° 
within the QCDF framework, which is based on the state- 
ment that the hadronic matrix elements of the operators 
in the effective weak Hamiltonian simplify in the heavy 
quark limit according to [l[ 



(MiM 2 \Qi\B) ~ F BAh (0)f M2 J duTl{u)<t> M M (1) 
+ /b/mi /m 2 / cLodvdu T? 7 (w, v, u)4>b M^Mi {v)4>m 2 (u). 



The factorization formula implies, on the one hand, that 
the theoretical prediction requires non-trivial hadronic 
input parameters, such as decay constants /, moments 
of light-cone distribution amplitudes <fi and form factors 
F, which encode all long-distance effects in the limit 
nit ~ * 00 ■ The power of the decomposition in ([I]) lies, 
on the other hand, in the fact that it provides the path 
to a systematic implementation of radiative corrections. 
The short-distance hard-scattering kernels T^' 11 are per- 
turbatively calculable and currently being worked out to 
next-to-next-to- leading order (NNLO) 0, la Q- 

The NNLO calculation is to date incomplete, but a 
subset of hard-scattering kernels, which specify the so- 



called topological tree amplitudes, has recently been de- 
termined to NNLO [4j, |(| . This allows us to present the 
first complete NNLO prediction within the QCDF frame- 
work for the decays B~ — » it~ ir / p~ p° , which are pure 
tree decays in the excellent approximation that small 
electroweak penguin amplitudes are neglected [28[ . 

As the considered decays are likely to be dominated 
by their standard model contribution, they may serve 
as benchmark channels for testing the various theoret- 
ical descriptions of the strong interaction dynamics in 
hadronic B decays. By normalizing the hadronic decay 
rates to their semileptonic counterparts B — > itiv/piv 
at maximum recoil, most of the theoretical uncertainties 
from hadronic input parameters and \V u b\ drop out and 
one obtains precision observables for testing the QCD 
dynamics of the topological tree amplitudes. We con- 
front the NNLO prediction in QCDF with experimen- 
tal data and find support for the factorization assump- 
tion. We also take a look at the other tree-dominated 
B — > irn/np/pp decay modes and conclude that the 
colour-suppressed tree amplitudes seem in general to be 
somewhat enhanced, which may hint at a smaller value 
of the first inverse moment of the B meson light-cone dis- 
tribution amplitude, Xb — 250 MeV. We conclude our 
analysis with a comment on the so-called B — > irir puzzle. 



II. TREE AMPLITUDES 

The decay amplitudes for hadronic B meson decays 
are conveniently parameterized by a set of topological 
amplitudes, which contain short-distance QCD and some 
electroweak effects. In the notation of |7j they read 
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with X p = V pb V* d and i„ = iG F /V2 m%f„F^ ^(0) and 
similarly for B~ — > p~ p° with f„ — > f p , F_ 
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TABLE I: List of input parameters (in units of GeV or di- 
mensionless) . Scale dependent quantities refer to p, = 1 GeV. 



0^(7777) — > ai(pp). Whereas the electroweak penguin am- 
plitudes E ,„ are currently known for B — > irir/irp/pp 

to NLO 0,[3, the tree amplitudes e_,2 have recently 
been determined for B — > 7777 to NNLO [1, From 
the projection properties of the leading-twist 77 and p 
wave functions, we find that the respective expressions 
for B — > ttp/plPl are identical (L refers to the longitu- 
dinal polarization). We in particular do not consider de- 
cays into transverse p mesons, which cannot be described 
model-independently as they do not factorize. 

We evaluate the tree amplitudes with 3-loop running 
coupling constant and NNLL Wilson coefficients @ of the 
operators in the weak effective Hamiltonian (we use the 
operator basis from (__]). The spectator scattering part 
(T/ ) receives contributions from two perturbative scales, 
p h ~ m b and p hc ~ (AQ C D"^f,) 1/2 , which gives rise to log- 
arithms lnmb/AQCD that we resum via renormalization 
group equations in SCET to LL approximation. Other 



scale dependent quantities are treated as described in [6j , 
except for the parameters of the B meson wave function, 
which we evolve with fixed order relations as their evolu- 
tion from their input scale does not induce parametrically 
large logarithms. 

We also include certain power corrections to the tree 
amplitudes that are related to subleading-twist wave 
functions of the light mesons. As these chirally enhanced 
contributions do not factorize, we use the model proposed 
in [l[ to estimate their size. 

Our theoretical input parameters are listed in Table [U 
We deduced our default values for the hadronic parame- 
ters from recent lattice and sum rule calculations (where 
available). In general the parameters related to the pion 
(fn, a£[Tl],G3, i^pj]) are better determined than the 
ones related to the rho meson (/pQ, <[H,[ll,^ p [ll]). 
While there exists a large number of calculations for the 
B meson decay constant Jb , less is known about the 
moments of the B meson wave function (Ar. en ?)p^|. 
Our value for As is based on a QCD sum rule calculation 
and on estimates from the operator product expansion, 
accounting for recent claims that higher dimensional op- 
erators lower the value of As (last paper of fill]). 

We estimate the size of higher order perturbative cor- 
rections by varying the factorization scales pu and phc 
independently within the ranges specified in Table [J 
On the other hand we evaluate the non-factorizable 
power corrections at a fixed scale po = 1.5 GeV. 
The latter introduce certain model parameters (pu , 4>h) 
and some additional hadronic parameters. We use 
(m u + m d )(2GeV) = 8 MeV, m b (m b ) = 4.2 GeV and 
/^(lGeV) = 165 MeV. 

This brings us to our NNLO prediction of the colour- 
allowed (ai) and colour-suppressed (c_) tree amplitudes. 
In the B — > 7777 ■/ ' plPl channels we obtain 
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where the uncertainties in the intermediate results stem 
from the variation of hadronic input parameters, higher 
order perturbative corrections and the considered model 
for power corrections, respectively, which have been 
added in quadrature for our final error estimate. 

We see, on the one hand, that the colour-allowed tree 
amplitudes ai can be computed precisely in the factor- 
ization framework. The colour-suppressed amplitudes c_ 
suffer, on the other hand, from substantial theoretical 



uncertainties. The problem is related to certain cancella- 
tions between various perturbative contributions, which 
make the real parts particularly sensitive to the specta- 
tor scattering mechanism which is proportional to the 
hadronic ratio /Mi/_/A_i ? ' BMl (0). Our poor knowledge 
of the B meson parameter A^ in particular translates 
into the uncertainties _o 049 anc ^ -0 043 f° r ^ ne rea ^ P ap ts 
of 02(7777) and Oi2(pLPL,), respectively. 
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III. BRANCHING RATIOS 

The branching ratios of B~ — > Tr~7r°/p~p° depend 
in addition on clectroweak penguin amplitudes, cf. |(3J). 
These amplitudes have not yet been determined to 
NNLO [29|], but their numerical values are rather small 
(l a 3/4Ewl ~ 0.01). As they are not CKM-enhanced in 
tree-dominated decays, it is consistent to treat these am- 
plitudes in the NLO approximation. The explicit NLO 
results can be found in [l|, @, Q (they are formulated in 
a different operator basis of the effective Hamiltonian). 
The CP-averaged branching ratios become 



10 6 Br(B- 
10 6 Br(B" 



7T-7T ) = 6.22 



+ 1.14 +2.03 +0.16 +0.43 
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01 n+8-5 
zl - u -7.3> 



(4) 



where the uncertainties in the intermediate results are 
due to CKM parameters, hadronic parameters, higher or- 
der perturbative corrections and non-factorizable power 
corrections, respectively. 

Our NNLO results are in good agreement with exper- 
imental data [19, ufll , 
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10 6 Br(S- 
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pIpD 
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= 22.5 



+0.41 
0.40' 
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i.e. the experimental values are reasonably well repro- 
duced by the central values of our NNLO prediction, 
which is based on the input parameters from Table|TJ One 
should keep in mind, however, that we could also have 
obtained similar numbers for the branching ratios with 
rather different values of the tree amplitudes, the form 
factors and \V u b\- As we discuss in the following section, 
a much stronger test of the factorization assumption can 
be obtained by considering the ratios of hadronic and dif- 
ferential semileptonic decay rates, where the dependence 
on the form factors and \V u b \ drops out to a large extent. 

We may also take a look at the other tree-dominated 
B — > mr /ftp/ pp decay modes. We emphasize that the 
NNLO calculation of these branching ratios is to date 
still incomplete, since the QCD penguin amplitudes have 
not yet been determined to NNLO (this is why we do not 
discuss CP asymmetries in this letter) . These modes also 
differ conceptually from B~ — ► Tr~ir°/p~p° in the sense 
that they receive contributions from weak annihilation, 
which constitutes another class of non-factorizable power 
corrections. We again use the model from |l[ to estimate 
their size. 

Our results for the CP-averaged branching ratios 
are shown in Table [TTJ Apart from some exceptions 
(7r + 7r~, 7r°7r°, 7r + p~, ir~ p + ) our default prediction (with 
central values) is again in reasonable agreement with 
the data. The agreement is, however, less pronounced 
than for the pure tree decays B~ — > tt^tt /p~p°. More- 
over, we point out that the colour-suppressed modes 



(7r°7r°, ir°p°, p°p°) are subject to sizeable theoretical un- 
certainties. This is partly related to the problem men- 
tioned at the end of the previous section (\b) and in 
addition to the fact that these modes are more likely to 
be affected by l/m^-corrections. 

In order to illustrate the correlation of the theoretical 
uncertainties, we show in Table |TT] the central values of 
some extreme scenarios (in the spirit of Q): 

In Scenario A we study the dependence on the weak 
phase 7 (we set 7 = 110°). Modes that show a strong 
dependence on this scenario (it + it~ ,ir°ir° , p° p°) are not 
particularly suited for our purposes, as we focus on test- 
ing the QCD dynamics of the topological tree amplitudes 
in this work. 

In Scenario B we pursue the question if the data are 
in accordance with a large colour-suppressed amplitude, 
which may be realized in the factorization framework by 
a very low value of Xb — 200 MeV (we moreover decrease 
the form factors to Ff ff (0) = 0.21 and A% p (0) = 0.27). 
This scenario shows a satisfactory description of the data, 
in particular the "problematic" modes 7r + 7r~, n + p~ and 
7r~/? + are - by construction - in much better agreement 
with the data. 

It is tempting to understand the large experimentally 
observed 7r°7r° branching ratio as an indication for size- 
able non-factorizable power corrections. It is hard to ad- 
dress this issue in a model-independent way. We would 
like to emphasize, however, that some observables are in- 
deed more likely to be affected by 1/ mb-corrections than 
others (cf. the column labelled "pow" in Table HT|) . We in 
particular expect the branching ratios of the tree decays 
B~ — > ir~n°/p~p a to be clean observables as they are 
free of weak annihilation contributions. 

In order to quantify this question we study the in- 
fluence of a large annihilation amplitude in Scenario C 
(within the BBNS model). It turns out that it is al- 
most impossible to enhance the 7r°7r° decay rate and 
to simultaneously decrease the 7r + 7r~ rate without fine- 
tuning the model parameters [301 ] . Moreover, the overall 
pattern of branching ratios and in particular the rates 
of the other colour-suppressed modes seem to disfavour 
a generic scenario with large annihilation contributions. 
This is illustrated in Scenario C, where we double the 
default value of the BBNS model for universal weak an- 
nihilation, i.e. we set pa = 1 and <f>A = 0. We conclude 
that we do not see any clear pattern of abnormally large 
power corrections in the data and prefer to be guided by 
clean observables rather than by the colour-suppressed 
and penguin-contaminated 7r°7r° branching ratio, which 
cannot be predicted precisely in the factorization frame- 
work anyway. We admit that our conclusion is a model- 
dependent statement, which is, however, supported by 
a light-cone sum rule analysis, which finds even smaller 
annihilation contributions than the BBNS model with 
default parameters [2l| . 

Finally, Scenario D is motivated by our analysis of the 
following section. It combines elements from Scenario A 
and B, but is based on a more moderate parameter choice 
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TABLE II: CP-averaged branching ratios (in units of 10 -6 ). Experimental values for B — » -kit /up are taken from fisl ] . whereas 
the ones for B — > pLpL have been inferred from [2pj ] . The different scenarios correspond to: large 7 (A), large colour-suppressed 
amplitude (B), large weak annihilation (C) and a combined scenario (D). Further details are given in the text. 



7 = 90°, X B = 250 MeV and Ff n (Q) = 0.23, which are 
within the ranges of our default parameters from Table U 
These values are inspired by a fit to a set of particularly 
clean observables that we discuss below. We refrain from 
presenting the details of our fit and prefer to simply il- 
lustrate the effects of such a combined scenario [31j ]. 



IV. PRECISION OBSERVABLES 

Our predictions for the branching ratios from Table [TT] 
typically have ~ 40% uncertainties, which are largely re- 
lated to an overall normalization from IV^I-F^^O) and 
|Kb|^ P (0). This particular source of uncertainties can 
be eliminated by normalizing the hadronic decay rates to 
the differential semileptonic rates at maximum recoil, 
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and similarly for B° — > p + i v\ with Ff 



Btt 



Aq P and 

— > m p . The situation is, however, different for the 
colour-suppressed modes (tt°tt° ,tt° p° , p° p°), which are 
rather dominated by the uncertainties from As and power 
corrections than by form factor uncertainties and \V ub \. 
We therefore do not consider these modes in this section. 

The BaBar collaboration has measured the semi- 
leptonic B — > irlv decay spectrum to high accuracy [22| . 
The data has been investigated in detail under different 
types of form factor parameterizations in [23| . This anal- 
ysis uses the HFAG average for the absolute branching 
ratio and finds \V ub \Ff "(0) = (9.1 ± 0.7) • 10" 4 , which 
is to be compared with our default value 10.3 • 10~ 4 and 
8.3 • 10~ 4 from Scenario B. The experimental value has 
been adopted in conjunction with our default value for 
\V ub \ to fix the form factor Ff 7r (0) = 0.23 in Scenario D. 



The analysis of the differential semileptonic B — > p£v 
decay spectrum is more complicated as three different 
form factors contribute in this case (which confine to 
I Vub I j4^ p (0) at maximum recoil). Recent measurements 
by BaBar. Belle and CLEO provide data in 3-4 q 2 - 
bins [24|], which does not yet allow to extrapolate the 
decay spectrum in a model-independent way. In a recent 
analysis the data has been combined with (quenched) lat- 
tice calculations of the form factors in the high q 2 region 
and light-cone sum rule predictions for q 2 = [25] . This 
analysis yields \V ub \A% p (0) = (5.5 ± 2.6) • 10~ 4 , which 
illustrates that the data is still premature. We therefore 
do not include this number in our analysis. 

Our predictions for the ratios 



Km 3 (MiM 2 ) 
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are shown in Table IIIII For the 7iy>modes we chose 
the normalization such that the dependence on the form 
factor multiplying the colour-allowed amplitude is most 
strongly eliminated. From Table IIIII it can be seen that 
the theoretical uncertainties have been reduced consider- 
ably to the level of ~ 15%. Moreover, correlations among 
different sources of theoretical uncertainties have been re- 
solved to a large extent. 

The first two ratios in Table|TTT]provide particular clean 
probes of the QCD dynamics of the topological tree am- 
plitudes 7, 26]. In the factorization framework we have 



ot\ + a 2 \ 



(8) 



where small electroweak penguin amplitudes have been 
suppressed. Our NNLO prediction for this ratio 

^(tt-tt ) = (0.70l° : 42 )GeV 2 (9) 

is in good agreement with experimental data 

7^-7r )| cxp - (0.8ll°:l 4 )GeV 2 , (10) 
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Theory 



CKM had 



/' 
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+0.04 
0.05 
+0.19 
0.14 
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18.7 



+0.08 
-0.07 
+0.16 
-0.18 
+0.15 
-0.17 
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+0.15 
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+0.31 
-0.31 
+0.03 
-0.04 
+0.10 
-0.07 
+22.7 
-10.8 



+0.17 

-0.09 

+0.18 

-0.12 

+0.03 

-0.06 

+0.15 

-0.19 

+0.03 

-0.06 

+0.16 

-0.21 

+0.13 

-0.07 

+0.13 

-0.13 

+0.02 

-0.02 

+0.14 

-0.07 

+7.0 

-11.0 



+0.03 
-0.03 
+0.03 
-0.05 
+0.02 
-0.02 
+0.05 
-0.06 
+0.02 
-0.02 
+0.06 
-0.07 
+0.03 
-0.03 
+0.00 
-0.00 
+0.00 
-0.00 
+0.03 
-0.03 
+2.6 
-2.3 



+0.11 
-0.08 
+ 0.11 
-0.08 
+0.16 
-0.09 
+0.23 
-0.17 
+0.18 
-0.10 
+0.28 
-0.22 
+0.08 
-0.08 
+0.06 
-0.14 
+0.02 
-0.00 
+0.08 
-0.10 
+ 10.2 
-10.4 



1.01 
1.35 
0.75 
2.44 
1.27 
2.61 
0.70 
2.06 
0.31 
0.90 
9.33 



1.16 
2.07 
0.97 
2.46 
1.01 
2.68 
0.89 
2.65 
0.39 
0.98 
8.32 



0.93 
1.79 
1.24 
2.99 
1.29 
3.22 
0.59 
2.49 
0.42 
0.57 
17.3 



1.07 
1.71 
0.86 
2.44 
1.13 
2.64 
0.79 
2.33 
0.35 
0.95 
8.13 



1 C7+U-32 
1 - J '-0.32 

o.8ot°:li 

2.431S-S 



0.89t°;i 4 4 



+0.69 
0.69 
+0.04 
0.04 



3.23 
0.33 



1 01 + 009 

1 - ul -0.09 



+0.43 
0.43 



3.33 



TABLE III: Ratios TZm 3 {M 1 M2) of hadronic and differential semileptonic decay rates as defined in (J7J) (in units of GeV 2 ) and 
ratios R(M\M2 / M3M4) of hadronic decay rates from (|11[) . The different scenarios A-D are described in the caption of Table ITT1 



which strongly supports the factorization assumption. 
It is, however, interesting that the central experimental 
value is in between our default prediction and the value 
0.95 GeV 2 from Scenario B, which may hint at a some- 
what larger value of the colour-suppressed amplitude and 
hence a lower value of the parameter Xb — 250MeV 
(which we adopt in Scenario D). Experimental data for 
the ratio lZ p (p~[p L ) could help to clarify this point. 

We recall that all other ratios from Table IIIII receive 
contributions from QCD penguin amplitudes that are 
not yet completely available to NNLO. Among these 
TZTr(n + p^), TZ p {p\ j p L ) and 1Z p (tt~p + ) are particularly 
suited to test the dynamics of the colour-allowed am- 
plitudes. Our prediction for TZ n (7T + p~) compares again 
well to the experimental value. 

The fourth colour-allowed ratio ^^(vr+Tr - ) is special, 
since the interference of the colour-allowed amplitude 
with the QCD penguin amplitude is not negligible in 
this case. This ratio is thus particularly sensitive to the 
choice of the weak phase 7. One should keep in mind, 
however, that the power corrections from weak annihila- 
tion represent another important source of uncertainties 
for this ratio. It is interesting to replace the BBNS model 
for weak annihilation by the light-cone sum rule predic- 
tion from [2ll ] . which strongly reduces the uncertainties 
from weak annihilation and hence enhances the sensi- 
tivity to 7 (we then find l-OSi^™^™ 2 !^ GeV 2 ). 
The current experimental value may then be considered 
as a hint at a large value 7 > 90°. A smaller value of 
7, on the other hand, may then imply the presence of an 
additional contribution to the QCD penguin amplitude 
or that power corrections, which are neither from weak 
annihilation nor from chirally enhanced wave functions, 
have been underestimated in our approach. The latter 
would be conceptually important, as it would increase 
the total uncertainty from power corrections in QCDF. 
We refrain, however, from drawing any conclusions con- 



cerning TZ^ItT^ 



and its implications for 7, as long 



as the penguin amplitudes have not been calculated to 
NNLO. 

In Table [TTT1 we also show some ratios of hadronic decay 
rates defined by 



R{M l M 2 /M^M i ) 



T(B -> M1M2) 
T(B' -> M3M4) 



(11) 



The ratio i?(p£ p\ I p\ P~l) yields complementary infor- 
mation on the tree amplitudes from the p-sector, where 
the contamination from the QCD penguin amplitudes is 
known to be less important [1, [27| • We consider the ex- 
perimental value for this ratio as another important evi- 
dence in favour of an enhanced colour-suppressed ampli- 
tude (Scenario B or D). 

The ratios R{p\ p^ / 'ir~ p + ) and R(tt + tt~ /7t+ p~) of 
colour-allowed modes can be predicted precisely in the 
factorization framework. Whereas the second ratio is in 
nice agreement with the data, the first one seems to some- 
what disfavour a scenario with a large weak phase 7. 

The last two ratios from Table HTT1 finally refer to what 
is known as the B — > 7T7r puzzle. Whereas the ratio 
R(ir~n° /ir + ir~) is by construction in Scenarios A, B and 
D in better agreement with experimental data than our 
default prediction, the ratio R(ir + ir~ /ir°ir°) illustrates 
what we mentioned at the beginning of this section, 
i.e. the bulk of theoretical uncertainties does not drop 
out in ratios that involve colour-suppressed modes. The 
uncertainties of our default prediction 



i?(7r + 7r~/7r 7r°) = 25.7 



+26.0 
18.7 



(12) 



are thus extremely large and the central value is in vast 
disagreement with the data. This ratio may be brought 
down by a factor of ~ 3 in Scenarios A, B and D, which 
may be seen as an independent evidence in favour of these 
scenarios. The fact, however, that these predictions still 



6 



suffer from ~ 60% uncertainties related mainly to the 
power corrections, af, Hhc and /b (in decreasing order 
of importance), shows that we cannot expect to predict 
this ratio precisely. We would like to add that there is no 
such puzzle in the irp/ pp channels, i.e. there is no general 
failure of QCDF to describe colour-suppressed modes. 

V. CONCLUSIONS 

We presented the NNLO QCDF prediction for the 
approximate tree decays B~ — > ir~it a / p~p° and up- 
dated the global analysis of the other tree-dominated 
B — > irir/irp/pp decay modes. Our analysis from Sec- 
tion [IV] showed that QCDF yields precise theoretical pre- 
dictions for particular ratios of decay rates. We find in 
general support for the factorization assumption and un- 
covered some hints for enhanced colour-suppressed am- 
plitudes, which translate in QCDF into a small value 
of the B meson parameter As- Theoretical progress 



from non-perturbative methods on the hadronic ratio 
/mi/b/AbF BMi (0) as well as experimental measure- 
ments of the semileptonic B — > ptv decay spectrum may 
shed further light on this issue. 
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